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Abstract. The chemical selectivity and great sensitivity of the Extended X-ray Absorption Spectroscopy 
technique allowed the determination, in the paramagnetic phase, of the structural distortions induced by 
doping in the spin-Peierls CuGeOa compound. The distorted environments were analyzed as a function of 
concentration, magnetic nature of impurity and the substitution site (Ni, Mn and Zn impurities on the Cu 
site. Si impurity on the Ge site). This has led to estimate the variation of the angles and pair distances, 
and hence to evaluate the magnetic coupling along the Cu chains in the vicinity of the impurities. The 
antiferromagnetic interaction between Cu first neighbors in the pure sample is found to be weakened around 
Ni, almost cancelled in the case of Mn doping, and even to change sign, producing a ferromagnetic coupling 
for Si doping. More generally, the structural distortions on a local scale are shown to be key parameters 
for the understanding of the magnetic properties of doped spin-Peierls compounds. 



PACS. 61.10.Ht X-ray absorption spectroscopy: EXAFS, NEXAFS, XANES, 
impurities in crystals; microstructure - 75.10.Pq Spin chain models 



1 Introduction 

The first inorganic spin-Peierls (SP) compound, CuGeOa,]!] 
has already been extensively studied for its original mag- 
netic properties. The possibility of substituting the cationic 
sites, Cu or Ge, by different magnetic or non-magnetic el- 
ements, has brought new information and rose new ques- 
tions. 

The structure of CuGeOa can be described as a stack- 
ing of CuOg octahedra and Ge04 tetrahedra along the 
c direction of the structure. Each O belongs to both a 
Ge tetrahedron and a Cu octahedron. The stacking of 
distorted octahedra (elongated along the direction of the 
apical 0(1) and with a planar distorted square of basal 
0(2) in the perpendicular direction) results in chains of 
magnetic Cu^+ ions located at the center of edge shar- 
ing squares of 0(2) running along the c axis (Fig. [J). At 
room temperature, this compound then consists of uni- 
form chains of spins 5=1/2, interacting by superexchange 
through 0(2) atoms leading to an antiferromagnetic (AF) 
Heisenberg-type coupling (Fig. [T]). A magnetoelastic in- 
stability produces a dimerization of the chains below a 
critical temperature Tsp of 14.25 K. This leads to a sin- 
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Fig. 1. Schematic structure of CuGeOa 



glet (non-magnetic) ground state separated from the first 
excited triplet state {S — 1) by an energy gap of 2 meV. 
Below Tsp, the magnetic susceptibility drops to zero with 



2 



V. Simonet et al.: Effect of structural distortions on the magnetism of doped spin-Peierls CuGeOs 



a temperature activated behavior, and diffraction experi- 
ments reveal new superlattice peaks . 

Doping, either with Zn, Ni, Mg or Mn on the Cu site [31 
|4l[5] or with Si on the Ge site [61(7] , results into a decrease 
of Tsp and, at a lower temperature, the Neel temperature 
Tjv, into a transition to an AF long-range ordering (LRO). 
Some of these impurities are magnetic, like Ni with spin 
1 and Mn with spin 5/2 [8], while the others, Zn, Mg 
and Si, are non magnetic. They also have very different 
ionic radii whose consequence on the magnetic properties 
will be discussed in Sec. IV and V. A systematic study 
of the magnetic properties of single crystals with various 
doping levels (CuGei_a:Sia;03 with < a; < 0.08, and 
Cui_j,MyGe03 with M = Zn, Mg, Ni and < y < 0.1), 
all synthesized and characterized in the same way, has 
allowed establishing a quasi-universal phase diagram for 
the spin-Peierls transition [5] . The variation of Tsp versus 
concentration was found identical for all impurities on the 
Cu site, including Mn, which has been studied since Ref. 
[S] (for < y < 0.02), and for Si impurities applying a 
scaling factor y 3x. This is illustrated in Fig. 2 of Ref. 
[S] and in the alternative representation of the present 
Fig. [2] where no scaling factor ?/ « 3a; is applied in view 
of the discussion of Sec. V.B. Surprisingly, doping outside 
the spin chains is thus about three times more efficient in 
destroying the spin-Peierls phase than doping within the 
chainsQ In contrast, T/v varies with the type of impurity: 
it is higher for Mn and Si-doping and slightly lower for 
Ni-doping, in comparison with the value, found identical, 
for the non-magnetic impurities Mg and Zn. A difference 
between Ge site and Cu site doping is observed concerning 
the T;v(y, x) transition line which is continuous in the first 
case and presents a kink between a dimerized AF phase 
at low doping and a uniform AF phase at larger doping 
in the second case [T^]. This behaviour is characteristic 
of a second-order and a first-order phase transition, re- 
spectively. Another particularity concerns the easy axis of 
magnetization in the AF phase: it is the a— axis for Ni 
doping and the c— axis for all other substitutions, which 
can be explained by the single-ion anisotropy of Ni^+ in a 
distorted octahedral environment [13| . 

These results must be considered in the light of the dif- 
ferent couplings along the Cu chains. The magnetic prop- 
erties are mainly driven by two kinds of AF couplings: 
A superexchange coupling Cu-0(2)-Cu between first Cu 



^ Note that, while for Ni, Zn and Mg, the results of Gre- 
nier et al. [5] are in rather good agreement with the work of 
other groups on crystals of different origin, a large spread of re- 
sults exists concerning Si-doping. The slope of Tsp{x) from the 
works of Katano et al. ^ and Hiroi et al. JXlj is smaller than 
that found by Grenier et al. while it is similar in the work of 
Wang et aZ. [11], leading to the scaling factors y = x, y = 1.5 x 
and y ^ 2 — 3 X respectively. However, we are more confident 
in the results of Grenier et al. concerning the Si-doping data 
since the weaker effect measured by other groups suggests that 
not all Si ions enter the CuGeOa structure and/or part of them 
segregates (as was already observed in highly doped samples). 
This leads to overestimate x and thus results in an apparent 
smaller scaling factor between y and x. 
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Fig. 2. (Temperature - concentration) magnetic phase diagram 
for CuGei_3;Si^03 and Cui-yMyGeOa, with M=Zn, Mg, Ni, 
and Mn. 



neighbors with an AF exchange constant Ji = 120 — 180 
K and a frustrating AF coupling between second neigh- 
bors Ja = 0.24 - 0.36 Ji via the Cu-0(2)-0(2)-Cu path 
}14pi5j (Fig. [T]). The inter-chain couplings are one order 
and two orders of magnitude smaller along the b and a 
axes, respectively [TS]. These inter-chain couplings induce 
a coherence between the phases of the distortion waves 
on different chains. In CuGeOa, the dimerization pattern 
is out of phase for adjacent chains in the b and a di- 
rections. The strength and nature of the magnetic inter- 
actions are strongly dependent on the Cu-0(2) and Ge- 
0(2) distances and even more drastically on the Cu-0(2)- 
Cu superexchange angle (denoted a). The Goodenough- 
Kanamori- Anderson rules [TB] state that the exchange is 
ferromagnetic (FM) for an angle of 90°. The deviation of 
the Cu-0(2)-Cu angle from this value (a=99° in CuGeOa) 
will drive the exchange towards antiferromagnetism |17j . 
However, Geertsma and Khomskii pointed out that cal- 
culations taking into account this sole effect still yield 
a slightly ferromagnetic exchange coupling T8]: The hy- 
bridization of Ge orbitals with the 2py ones of the 0(2), 
which destroys the equivalence of 2px and 2py 0(2) or- 
bitals, must indeed be taken into account to obtain, for 
Ji, an AF exchange with a reasonable value of 135 K. 

The aim of the present work is to understand the effect 
of doping the Cu chains, with magnetic and non-magnetic 
impurities, and out of the chains in the case of Si-doping, 
by measuring the induced structural distortions of doped 
CuGeOa. These distortions should strongly influence the 
magnetic properties, in particular via the variation of the 
superexchange angle a and thus of the related magnetic 
coupHng Ji. The EXAFS (Extended X-ray Absorption 
Fine Structure) spectroscopy technique was chosen for this 
purpose as a local order probe, with chemical sensitivity, 
and for its possibility to measure very diluted systems. For 
the first time with this technique, the influence of Ni, Mn, 
Zn and Si-doping was investigated in the paramagnetic 
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Fig. 3. Comparison of the measured x(^) functions (a) for 
pure CuGeOs at the Cu K-edge, and for 6% Ni, 1.5% Mn and 
4% Zn doped samples, at the impurity A"-edges and (b) for 
pure CuGeOa at the Ge K-edge and for 2% Si-doped sample 
at the Si K-edge. 



phase (T > Tsp) but not that of Mg-doping because of 
the too low energy absorption of the Mg atom. 

The next section is devoted to the experimental details 
and to a brief recall of the EXAFS analysis methodology. 
The EXAFS results are reported in section III for Cu site 
and Ge site doping. The structural distortions, as deduced 
from the EXAFS refined parameters, are analyzed in sec- 
tion IV. Their influence is discussed in section V, in re- 
lation with the magnetic properties, and a conclusion is 
given in Sec. VI. 



2 Experiment 

The EXAFS measurements at the Cu and Ge absorption 
i^-edges, in transmission detection mode, and at the Ni, 
Mn and Zn absorption K-edges, in fluorescence detection 
mode, were carried out at LURE (Orsay) on the D42 and 
D44 beamlines of the DCI storage ring. The experiments 
were performed at 78 K, in the paramagnetic phase, above 
the spin-Peierls transition. The incident energy was varied 
by 2 eV steps using a Si(lll) single crystal monochroma- 
tor. The EXAFS oscillations were obtained in transmis- 



sion detection mode by measuring the X-ray beam inten- 
sities with ionization chambers, before (/q) and after (Ji) 
the sample, in an energy range close to the absorption 
K-edge of Cu (Ag"=8979 eV) and Ge (A^'==11103 eV). 
The samples studied in the transmission detection mode 
were pellets of homogeneous thickness made of a mixture 
of about 15 mg of single crystal ground into very fine 
powder and 35 mg of cellulose. In the fluorescence detec- 
tion mode, the powder obtained from ground doped crys- 
tals was spread on an adhesive tape. The EXAFS spectra 
were measured in the fluorescence detection mode with a 
niultielements detector around the impurities A'-edges in 
the 2% and the 6% Ni-doped samples (absorption K-edge 
A^^ = 8333 eV), 1.5% Mn-doped sample (A^" = 6539 
eV), and the 4% Zn-doped sample (Af" = 9659 eV). 

The Si K-edge EXAFS experiment was performed on 
the SA32 beamline of SUPERACO using a 2 InSb crystals 
monochromator. A thin slice of 2% Si-doped single crys- 
tal (absorption K-edge A|? = 1839 eV) was cleaved and 
stuck with Ag glue on an indium covered sample-holder. 
The absorption signal was recorded with a monoelcment 
detector in the fluorescence detection mode. 

The EXAFS oscillations were analyzed using the FEFF 
package [19^. The EXAFS oscillations, xif^)i were extracted 
from the experimentally measured absorption coefficient 
fi{E) using the AUTOBK program with k, the wavenum- 
ber of ejected photoelectron, equal to: 



2m{E- Eo)/h^ 



(1) 



where Eq is the absorption energy for the if -edge and m 
the electron mass. 

For an absorbing central atom surrounded by N iden- 
tical atoms, when considering only single scattering paths, 
EXAFS oscillations created by the backscattering of the 
photoelectron from the neighboring atoms are given by: 

xik)-~S'o{^ |F(fc,7r)|e-2'=^'^^-2^AW 

sin [2fcr -1- 2(5(fc) -f <?>(fc)] (2) 

where r is the mean pair distance, is the associated 
Debye- Waller factor which takes into account the struc- 
tural and dynamical disorders, and Sq is the amplitude 
reduction factor accounting for many-body effects within 
the absorbing central atom. The other parameters in Eq. 
([2|) are the phase ^(fc) and the modulus of the complex 
function F{k, tt), which both depend on the nature of the 
backscattering atom, the mean free path X{k) and the 
absorbing atom phase shift d{k). They are derived from 
the theoretical ab-initio multiple scattering calculations 
by the FEFF program using crystallographic data. In or- 
der to describe a real environment, the surrounding atoms 
are gathered into atomic shells (atoms at approximately 
the same distance from the absorber), themselves subdi- 
vided into several subshells, for each kind of chemically 
and crystallographically equivalent neighbors. 

Representative measured EXAFS oscillations x(fc) are 
shown in Fig. [31 Note that the spectra of the doped sam- 
ples, recorded at the impurity K-edge, are more noisy than 
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Fig. 4. Structure of CuGeOa . The Cu and Ge environments are 
materialized: starting from the dashed Cu and Ge absorbers, 
the atoms labelled "A","B", ... "F" (for Cu) and "a", "b", ... 
"g" (for Ge) are their neighbors going from the closest to the 
furthest ones that were accessible in the EXAFS experiment. 



those recorded at the Cu and Ge K-edges because of the 
small amount of absorbers and of the fluorescence detec- 
tion mode. In consequence, the EXAFS analysis could not 
be performed up to such large k values at the impurity 
i^-edges compared to the Cu and Ge ones. For the Si K- 
edge spectrum, an additional source of noise was the low 
fluorescence yield and the small photon flux at this low en- 
ergy. The complex Fourier transforms (denoted FT) of the 
weighted x(fc) functions were then calculated using a 
Hanning window. The amplitude of this FT can be associ- 
ated with a pseudo-radial atomic distribution around the 
absorber in the r-space. Lastly, the EXAFS oscillations of 
the atomic shells of interest have been isolated by a subse- 
quent inverse Fourier-transform (denoted IFT) performed 
in a restricted appropriate r-range and leading to a Fourier 
filtered fc"^x(/c) function. The FT and IFT were simulated 
with the FEFFIT program by summing the contributions 
given in Eq. ^ for all subshells obtained from the struc- 
tural model. The multiple scattering paths, whose contri- 
bution was checked to be weak, were not included in the 
calculation. In the fit, the pair distances, the Debye- Waller 
factors and the energy shift AEq were allowed to vary0. 
For the analysis of the EXAFS oscillations recorded at the 



^ AEo accounts for inaccuracies in the FEFF determination 
of the absorption threshold energy Eq. Different values of AEq 
were refined for the O, Cu and Ge neighbors. These extra shifts 
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Fig. 5. EXAFS results for the pure CuGeOs sample. Ampli- 
tude of the FT of k^x{k) (left) and Fourier filtered k'^xik) 
(right) at the Cu (top) and Ge (bottom) K-edges. The mea- 
sured spectra are represented by full circles, the fit by a thick 
solid line. The contributions of each kind of neighbors are also 
shown on the FT plots (left). 

impurity absorption edges, the FEFF files for each scat- 
tering path were recalculated with the impurity replacing 
Cu or Ge absorber within the same structure. 



3 Results 

3.1 Pure sample 

The structure of CuGeOs in the paramagnetic state was 
initially refined by VoUenkle et al. [20 and later by Braden 
et al. JTTU^ (space group Pbmm, a = 4.796 A, 6 = 
8.466 ATc = 2.940 A at 295 K), from Rietvelt refinement 
of X-ray and neutron powder diffractograms 0- The crys- 
tallographic data of Braden et al. |17] were used as input 

are meant to correct errors introduced by FEFF in the calcu- 
lation of the scattering phase shifts. 

The structure was revised by M. Hidaka et al., J. Phys.: 
Condens. Matter 9, 809 (1997). They found a different space 
group P2i2i2i and a larger unit cell 2a x 6 x 4c than Braden et 
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of the FEFF program in our EXAFS analysis. These struc- 
tural parameters yield, for the atomic distribution around 
Cu atoms, a first shell with the 4 square planar 0(2) at 
1.942 A and a second shell including the 2 apical 0(1) 
at 2.78 A, the 2 chain neighbors Cu at 2.95 A and 4 Ge 
at 3.30 A (top of Fig.|4]). The atomic distribution around 
Ge atoms is made of a first shell with the 4 O of the tetra- 
hedron (2 0(2) at 1.737 A and 2 0(1) at 1.778 A) and 
a second shell including 2 Ge at 2.946 A and 4 Cu at 
3.30 A (cf. bottom of Fig. H]). This sheU descript ion IS 
summarized in Table I. 

The EXAFS fitting procedure was tested on the pure 
CuGeOs sample. The first two atomic shells of the ra- 
dial distributions at both Cu and Ge K-edges were simul- 
taneously fitted. This implies that the distance and the 
Debye- Waller factor corresponding to each Ge-Cu paths 
were hold identical for the two sets of data during the re- 
finement. The k range used for the FT and for the fitting 
procedure was [3.5, 15 A-i] while the r range used for 
the IFT was [1, 3.4 A] at both edges. The fits are in good 
agreement with the experimental spectra, as shown in Fig. 
[5l where the contribution of each kind of neighbors is rep- 
resented. Note that, at the Cu K-edge, the signal from the 
two apical oxygens 0(1) is already very weak and the fur- 
ther 0(1) neighbors were not considered because yielding 
a negligibly small contribution. However, the presence of 
the next 4 Ge atoms at about 3.786 A was found necessary 
to fit correctly the second shell. The resulting structural 
parameters obtained at both edges are listed in Table I 
and the comparison with their initial values shows the 
consistency of the fit. After this first EXAFS analysis, the 
same procedure was applied to the doped samples. 
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3.2 Substitution on the Cu site 

The Ni and Mn i^-edges EXAFS oscillations were ana- 
lyzed up to fc = 12 A~^ (see Fig. |3|). In the case of Zn, an 
additional difficulty arose from the mixing of Cu K-edge 
and Zn K-edge EXAFS oscillations due to the closeness 
of the absorption energies of both elements: 8979 eV for 
Cu and 9659 eV for ZnQ Therefore, the Zn K-edge spec- 
trum, after correction from the contamination by the Cu 
K-edge oscillations, could be analyzed in a more limited 
k range, up to 9 A"^ (cf. Fig.O. 

First, the structural results derived at the Cu K-edge 
in the doped samples were found identical to those of the 
pure sample, even for doping levels as high as 6% (Ni). 
The distortions induced by the impurities were then de- 
duced from the comparison of the structural parameters 



al. |17| . who nonetheless reconfirmed the initial structure one 
year later \21\ . 

* To correct the Zn K-edge spectrum, the Cu isT-edge oscil- 
lations x{E) of the pure sample, recorded up to 9998 eV, were 
converted into x{k) using Eq.([T| with the absorption energy 
Eq of Zn. The Cu EXAFS contribution, whose amplitude at 
the Zn edge was obtained from the ratio of both absorption 
jumps recorded in a unique scan in the doped sample, was 
then subtracted from the Zn spectrum. 



Fig. 6. EXAFS results for the 2% and 6% Ni-doped CuGeOs. 
Amplitude and imaginary part of the FT of k^x{k) (right) 
and Fourier filtered k'^x(k) (left). Comparison of the measured 
spectra at the Ni if-edge for both doping levels and at the Cu 
K-eAge for 6% Ni sample (top). Comparison of the measured 
spectra (circles) and of the fits (solid lines) at the Ni K-edge 
for 2% doping (middle) and 6% doping (bottom). 



obtained at the impurity and at the Cu i^-edge in the 
same sample and using the same k and r-ranges of anal- 
ysis for the EXAFS oscillations. Note that we did not try 
to take into account impurity neighbors in our structural 
model (although there can be some for high doping lev- 
els) because it is most probably impossible to distinguish 
them from Cu neighbors in the EXAFS analysis. 

The FT of fc'^x(fc) and Fourier filtered k^xi^) ^re shown 
in Fig. [H] for the 2% and 6% Ni-doped compounds and in 
Fig.Elfor the 1.5% Mn and 4% Zn-doped compounds. The 
spectrum recorded at the impurity isT-edge in each com- 
pound is compared to its fit and to the reference spectrum 
at the Cu K-edge. The structural parameters extracted 
from the fits are reported in Table II. The fits in k and 
r spaces are of rather good quality. For all doping values, 
the amplitude of the FT, for the first O shell, is smaller at 
the impurity if -edge than at the Cu one, which is related 
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Table 1. Structural parameters of the atomic radial distribution up to 3.95 Acentered on Cu and Ge atoms, from the crystal- 
lographic data of Braden et al. ITTI and from the EXAFS simulations of the pure CuGeOs sample at both Cu and Ge if-edges. 
The atom labels are those of Fig. |4] Crystallographic data of Braden et al. [17] imply two slightly different Ge-0 distances for 
the two pairs of oxygens within the tetrahedron. The parameters reported with no error bar were held fixed during the fit. The 
error bars only account for the statistical uncertainties and are therefore probably underestimated. The pair distances r and 
Debye- Waller factors are in A and A^ respectively. The quality of the EXAFS fit is given by the 7l-factor: 71/ =0.002 at the 
Cu Tf-edge and Rf =0.008 at the Ge if-edge. The refined Sl values are 0.927 at the Cu Tf-edge and 0.953 at the Ge K-eAge. 
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Zn [Ar/r = 4.5%) and Mn {Ar/r = 7.3%). Interestingly, 
the structural distortion varies also for a same impurity as 
a function of its doping level, as shown by the increase of 
the Debye- Waller factor and first Ni-0(2) distance from 
2% to 6% Ni. For Ni and Zn impurities, there is an ambigu- 
ity concerning the position of the next nearest neighbors, 
the apical oxygens 0(1). They are found either closer or 
further away from the impurity than from Cu (see param- 
eters of the equally good EXAFS fits in Table II). At the 
Mn K-edge, there is no such ambiguity since the EXAFS 
spectra is 10 times better fitted when the 0(1) atoms are 
closer from the impurity than further away, as compared 
to the Cu environment. Lastly, there is a clear variation 
of the distances to the next Cu and Ge neighbors with 
respect to the Cu environment (1.6% and 1% increase, 
respectively) in the case of the Mn impurity only. 



3.3 Substitution on the Ge site 



Fig. 7. EXAFS results for 1.5% Mn (top) and 4% Zn (bot- 
tom) doped CuGeOs : amplitude and imaginary parts of the 
FT of k^x{k) (right) and Fourier filtered k^x{k) (left) of the 
measured data (full circles) and their fits (solid lines) at the im- 
purity K-eAge. The Cu iC'-edge spectra measured on the same 
sample are also shown (dashed lines), for comparison. 



to an increase of the Debye- Waller factor, in agreement 
with a doping-induced structural disorder. An increased 
Debye- Waller factor is also observed for the second shell 
in the case of the 6% Ni and 1.5% Mn-doped compound. 
Note that the 4% Zn-doped compound shells are less well- 
defined because of the smaller k range of integration. For 
all compounds, the first shell is shifted towards larger r 
values. The distance from the absorber to the 0(2) first 
neighbors increases more and more when replacing Cu by 
Ni {Ar/r = 1.4% and 2.2% for 2 and 6% doping resp.). 



The EXAFS oscillations, recorded at the Si K-edge in a 
2% Si-doped sample, are limited to 8.7 A~^ and allow get- 
ting information only about the first O tetrahedral shell. 
Note that, due to the small range of fitting, the Debye- 
Waller factor could not be accurately refined and was fixed 
to 0.002 A^, i.e. to the value found in the pure sample at 
the Ge K-edge. This may be one reason for the small dis- 
crepancies (i?-factor value of 0.011) observed between the 
fit and the experimental data (Fig. [8|). The results, ob- 
tained using 5q=0.81 in the fitting procedure, indicate a 
strong contraction of the Si-0 distance of about 8.4% with 
respect to the Ge-0 distance. In addition, a more symmet- 
ric environment is found around the Si than around the 
Ge atoms, that can be described by a small regular O 
tetrahedron with distance from center to corner of 1.61(2) 
A (instead of 1.74 and 1.78 A for the two pairs of oxy- 
gen forming the distorted tetrahedron surrounding the Ge 
atoms). 
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Table 2. Structural parameters obtained from the simulations of the EXAFS spectra at the impurities iC'-edges of the Cu- 
site doped CuGeOa samples: i?- factors (quality of the fit), relative distance variations Ar jr and Debye- Waller factors . The 
relative distance variations induced by impurities are calculated with respect to the results obtained for the Cu A'-edge spectrum 
of the same sample analyzed using the same k and r ranges (see Table [!}. For Ni and Zn impurities, two roughly equivalent 
solutions of the EXAFS fits are reported with a different distance from the impurity to the apical oxygens 0(1), while only one 
solution is relevant for Mn. The results given in bold characters correspond to the solutions validated by the analysis presented 
in Sec. 14.11 Note that the parameters presented with no error bar were held fixed during the fit and equal to the Cu is'-edge 
refined value in the same sample. The Sq factors at all the impurities Tf-edges were fixed to the value of 0.927 obtained at the 
Cu /f-edge. 





2% 


Ni 


6% 


Ni 


4% 


Zn 


1.5% Mn 




0.0021 


0.0023 


0.0025 


0.0058 


0.0015 


0.0030 


0.0016 




0(2) 


1.4(7) 


1.3(6) 


2.2(7) 


2.2(7) 


4.5(1.8) 


3.9(1.3) 


7.3(6) 




0(1) 


-17(4) 


4(4) 


-15(2) 


4(4) 


-15(3) 


9(2) 


-13(1) 


^(%) 


Cu 


0.6(6) 


0.7(7) 


-0.3(6) 


0.1(8) 


-0.9(8) 


-0.7(1.0) 


1.6(6) 




Go 


0.0(4) 


0.15(70) 


0.1(5) 


0.4(7) 


0.3(1.6) 


1.3(1.2) 


1.0(6) 




Ge 


0.0 


0.15 


0.1 


0.4 


0.3 


1.3 


1.0 




0(2) 


0.0038(4) 


0.0039(4) 


0.0042(4) 


0.0042(5) 


0.0035(11) 


0.0044(6) 


0.0045(3) 




0(1) 


0.025(12) 


0.017(13) 


0.0119(4) 


0.013 


0.017(13) 


0.01 


0.0088(20) 


(A^) 


Cu 


0.0034(8) 


0.0032(9) 


0.0038(6) 


0.0035(8) 


0.002 


0.002 


0.0047(10) 




Gc 


0.0032(4) 


0.0034(5) 


0.0045(4) 


0.0049(6) 


0.0045 


0.0045 


0.0048(4) 




Go 


0.021(11) 


0.025(19) 


0.012(3) 


0.014(6) 


0.014 


0.014 


0.026(16) 



tion of the first shell around the impurities with the pair 
distances evaluated from these ionic radii (tabulated in 
Table I of Ref . [22j ) . This qualitative analysis allows check- 
ing the validity of the EXAFS analysis and lifting some 
remaining ambiguities. 

In practice, the A/-0 pair distance (M denoting the 
absorbing ion on the Cu or Ge site) is evaluated by sum- 
ming the ionic radii of M and of its first neighbors 
for comparison with the distances obtained from EXAFS. 
The results are summarized in Table IIL The simplest 
case concerns the Ge site in which Ge'*"'" and Si^+ ions are 
surrounded by a slightly distorted and a regular tetrahe- 
dron, respectively. The calculations of the Ge-0 and Si-0 
distances determined from EXAFS are in excellent agree- 
ment with the distances derived from the ionic radii [^ . 

The case of the Cu site is more complex. The Cu^+ 
is surrounded by an elongated octahedron made of a dis- 
torted basal square of 4 0(2) and 2 apical 0(1) at a much 
larger distance. The calculated Cu-0 distance from the 
ionic radii, in the hypothesis of square planar configura- 
tion, yields a very good agreement with the measured Cu- 
0(2) distance indicating that the orbital configuration of 
the Cu is mainly sensitive to the 4 closer 0(2). Accord- 
ing to the EXAFS results, there are two indistinguish- 
able solutions for the geometry of the oxygen octahedra 
around Ni^+ and Zn^"*" impurities on the Cu site. The api- 
cal 0(1) would get closer to the impurity or further away 
as compared to the Cu environment, corresponding to an 
environment closer to a regular octahedron or closer to 
a square planar environment, respectively. Note however 
that a square planar configuration exists only for Cu and 
Ni. Concerning the Ni, the only way to reproduce the ob- 
served increase of the Ni-0(2) distance, as compared to 
the Cu one, is an orbital configuration compatible with 
an O octahedron around the Ni and excludes a square 




I 2% Si I 
• Si K-edge 




pure CuGeO , Ge /<-edge 



15 



2 4 6 8 10 12 140.0 0.5 1.0 1.5 2.0 2.5 
k(K') r{A) 



:15 



Fig. 8. EXAFS results for 2% Si-doped CuGeOs. Amplitude 
and imaginary part of the FT of fc'^x(^) (right) and Fourier 
filtered k'^xi^) (left) at the Si K-edge of the doped sample 
(full circles) compared to the fit (solid lines) and to the results 
obtained at the Ge A'-edge on the pure sample (dashed lines) . 



4 Analysis 

4.1 Comparison with ionic radii 

The ionic radius for each chemical species, which depends 
mainly on its oxydation and spin states, on its coordi- 
nation number (the radius increases with the number of 
neighbors) and on the polyhedral description of its envi- 
ronment, can be roughly estimated as a function of these 
sole criteria and independently of the detailed structure 
[22]. It is therefore interesting to compare the results of 
the present experiment concerning the structural distor- 
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Table 3. Distances to the first O neighbors, in A, calculated from the ionic radii from Ref. [22]) in octahedral and in 
square planar environments (when possible) for the Cu site, and in tetrahedral environment for the Ge site. The distances in 
bold are compatible with those derived from EXAFS and reported below {'') with the corresponding elongation ratio of the 
octahedron. The 2% and 6% Ni-doping results are reported. The last line of both tables gives the superexchange a angle in 
degrees, estimated from these EXAFS results (see Sec. 14. 2p . 





Cu 


Ni 


Mg 


Zn 


Mn 


Ge 


Si 


planar square" 


1.92 


1.84 












octahedron" 


2.08 


2.04 


2.07 


2.09 


2.18 






tetrahedron" 












1.74 


1.61 


0(2)" 
0(1)" 

r{M/Cu-0{l)) 
r(M/Cu-0(2)) 


1.942 


1.96/1.98 




2.03 


2.08 


1.737 


1.61 


2.78 
1.43 


2.26/2.33 
1.15/1.18 




2.36 
1.16 


2.38 
1.14 


1.778 


1.61 


a angle 


98.9(1) 


97.9(3) 




95.8(3) 


95.85(15) 


98.9(1) 


94(1) 



planar description. This conclusion also holds for Zn, for 
which only an octahedral environment is expected to ex- 
ist. For Mn2+, the EXAFS analysis was unambiguous. A 
Mn-0 distance close to the one derived from the EXAFS 
analysis is indeed found when the ionic radius of Mn^+ is 
evaluated assuming an octahedral environment and a high 
spin state (i.e. S=5/2), as inferred from the fine structure 
splitting of the ESR spectra at low temperature [8] . Note 
that the M-0 distances for Ni, Zn and Mn, are slightly 
smaller than the one derived from the ionic radii because 
the octahedron is distorted. Some characteristic features 
of elongated octahedral environments are also observed in 
the X-ray Absorption Near Edge Structure spectra mea- 
sured at the Cu K-edge in the pure compound and at 
the Cu-site impurities K-edge in the doped compounds, 
in agreement with the present analysis. These considera- 
tions allow to validate the EXAFS solutions for the M- 
O distances and to reproduce qualitatively their increase 
when substituting Cu by Ni, Zn and Mn, in this order. For 
the non-magnetic Mg impurity, extensively studied but for 
which no EXAFS measurements could be done, the follow- 
ing prediction can be made based on the same ionic radii 
analysis: the distortions around Mg are expected to be in- 
termediate between those induced by Ni and Zn. Note that 
the distortion reduction of the oxygen octahedron, when 
replacing Cu by the other impurities, could be expected 
since the axial elongation in the Cu^"*" case is related to 
a strong Jahn- Teller effect, usual for this ion. Its substi- 
tution with non-magnetic or non Jahn- Teller 3d ions, like 
Ni^"*" or high spin Mn^+, should reduce this effect. The 
remaining octahedral distortion is then only induced by 
the structural environment. 



To summarize, i- the Cu environment is closer to a 
square planar one {rcu-O(i) /''^Cu-0(2) — 1-4) while sub- 
stituting this site yields an environment closer to a reg- 
ular octahedron (rM-o(i)/''A/-o(2) —1.14-1.18 for M — 
Ni, Zn, Mn) with the first distance increasing with re- 
spect to Cu in the following order: M = Ni, Mg, Zn, Mn, 
it- Si-doping on the Ge site induces a contraction of its 
surrounding oxygen tetrahedron, as strong as the pushing 
away of the 0(2) first neighbors around the Mn impurity. 



4.2 Superexchange angle calculation 

The distance variations in the impurities environment de- 
termined from the present experiment can now be used 
to calculate a structural parameter relevant for the mag- 
netic properties: the superexchange angle a, connecting 2 
adjacent Cu sites through the 0(2) atom. In the case of 
doping on the Cu site, this angle a = M— 0(2)— Cu will 
be modified within one chain on both sides of the impurity 
M, while in the case of doping on the Ge site, this angle 
a = Cu— 0(2)— Cu involves the 0(2) connected to the Si 
impurity and will thus be altered on two adjacent chains 
along the b direction (see Fig. [4|). 

Concerning the Cu-site impurities M = Zn, Ni and 
Mn, two hypothesis for the direction of displacement of 
their nearest neighbors can be made in order to estimate 
a from the measured r]^[_o{2) ^^nd rM-Cu distances. Note 
that the distance tm-Cu is altered only in the case of Mn 
doping. In the first hypothesis, the 0(2) are pushed away 
from the impurity in the same direction than the initial 
Cu-0(2) one (cf. Fig.[9ja), arrow (1)). This also produces 
a slight modification of the distance between the 0(2) and 
the first Cu neighbor of the impurity. In the second hy- 
pothesis, this Cu-0(2) distance is held identical to that 
in the pure sample, provided a small modification of the 
M-0(2) direction is made (cf. Fig.IlJa), arrow (2)). This 
assumption is suggested by the rigidity of the Cu-0(2) 
bond whose length is found unchanged in the otherwise 
largely distorted CuSiOa compound p3[|24] . The calcu- 
lated a values with both hypothesis are 98-97.7° for 2% 
Ni, 97.6-97.2° for 6% Ni, 96.1-95.5° for 4% Zn, and 95.7- 
96° for 1.5% Mn respectively, yielding the average values 
listed in Table III. The angle decreases for all impurities 
on the Cu site, with respect to their value in the pure 
compound (99°), but this effect is smaller in Ni-doped 
compounds than in the Zn and Mn-doped compounds. 

In the case of substitution by Si on the Ge site, a com- 
plication arises from the fact that the whole geometry is 
not planar i.e. the fi angle between the Cu-0(2) ribbons 
and the Ge-0(2) segment is equal to 159.5° (cf. Fig.IlJb)). 
A simple assumption for the Si-induced distortion is a con- 
traction of the O tetrahedron around the impurity without 
change of bond direction with respect to the Ge environ- 
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(1) 




Fig. 9. Schematic view of the local deformation, due to A//Cu 
(a) and Si/Ge (b) substitution, assumed for the calculation of 
the local Cu(or M)-0(2)-Cu angle a (see text for the different 
scenarii) . 

ment. In this case, the 0(2) atom is pulled towards the 
Si atom by 0.13 A (arrow in Fig. EJb)). This leads to an 
increase of the /? angle (161.7°) and a decrease of the a 
angle down to 94.1°. Note that a different variation of (3 
towards larger (resp. smaller) values would increase (resp. 
decrease) the value of a (estimated variation within ± 
1°). A possible contraction of the Cu-Cu distance was not 
taken into account in this simple estimation of a, while 
it is suggested from the study of Si-doped polycrystals 
by X-ray diffraction (variation of the c parameter) [25 . 
This would decrease even more the Cu-0(2)-Cu angle. In 
conclusion, the reduction of the a angle should be even 
stronger for Si than for Zn and Mn-doping. 



5 Discussion 

5.1 State of the art 

The next step is to relate the EXAFS results, and espe- 
cially the above determination of the superexchange angle, 
to the magnetic properties of doped CuGeOa. Let us recall 
that there is a universal decrease of Tsp versus the doping 
concentration with y Ki 2>x m (Cui_j^My)(Gei_2,Sia;)03, 
whereas the lower temperature onset of an AF LRO oc- 
curs at different Tjv values depending on the impurity 5 . 
Two observations have stimulated a lot of theoretical work 
[26 l l27 p 28 p29U30p 31j: the coexistence of both AF and SP 
LRO for small Si, Zn, Ni and Mn doping, as shown by 
neutron scattering studies [32ll33l l4l[34] . and the apparent 
absence of critical concentration of impurities for the on- 
set of this AF LRO [3S]. From the different models that 
have been proposed, a generally admitted description of 



the magnetic excitation of these dimerized chains emerges 
in terms of topological defects connecting regions of oppo- 
site dimerization parity. These defects, that can be intro- 
duced through doping, are then constituted by a soliton, 
pinned to the Cu substituent, followed on the same chain 
by an antisoliton which restores the interchain phasing. 
Although mobile in the isolated chains, the antisoliton 
remains confined close to the impurity site when inter- 
chain couplings are accounted for in order to minimize 
the size of the defective segment which presents an out- 
of-phase dimerization with respect to the 3D pattern |361 
28,30,31^. Such a description of the solitonic defects, in 
terms of soliton-antisoliton pairs carrying an uncompen- 
sated spin 1/2, leads to large AF fluctuations coexisting 
with underlying reduced dimerization. In pure samples, 
the spin-Peierls phase is established when the critical size 
and phase coherence of the dimerized domains is reached, 
allowing the system to collapse into a long-range ordered 
state. This average size of coherent dimerized domains 
is limited by thermally activated and/or doping- induced 
solitonic defects. This qualitatively explains the decrease 
of Tgp with increasing impurity concentration. On the 
other hand, in the presence of inter-chain coupling, the on- 
set of an AF LRO occurs when the AF correlated regions 
associated with these doping-induced soliton-antisoliton 
pairs overlap. 

The magnetic behavior under doping is thus qualita- 
tively explained for non-magnetic impurities within the 
Cu chains. They simply cut one dimer, thus releasing one 
spin 1/2 solitonic defect confined close to the impurity. 
However, this is less obvious in the case of magnetic im- 
purities on the Cu site, which may still be magnetically 
coupled to the neighbor Cu spins, and even less for Si im- 
purities that are introduced outside the Cu chains. The na- 
ture of the released soliton-antisoliton pair through doping 
is a key feature in the understanding of the establishment 
of the AF LRO. Its nature will crucially depend on the 
strength of the local magnetic interaction. It is therefore 
necessary to relate the structural parameters determined 
in the vicinity of the impurity in the previous sections with 
the value of the Cu-M exchange interaction. The variation 
of the in-chain nearest neighbor constant, Ji, with the 
Cu-0(2)-Cu superexchange angle a was inferred from the 
pressure dependence of the susceptibility in CuGeOa and 
found to be equal to: 5Jx/{Ji5a) w5.8% [37l[38| . Various 
calculations based on the analysis of microscopic magne- 
toelastic coupling 39j and microscopic structural mod- 
els [17] yielded different estimations: 6Ji/{JiSa) ~ 10% 
and 22%, respectively. A value of 6Ji/{Ji5a) ranging be- 
tween 10 and 20% was also found in microscopic calcula- 
tions of the exchange constants carried out using different 
approaches: perturbation theory, exact diagonalization of 
CU2O2 clusters, and band calculations [38]. The general 
result of these studies is that a reaches its critical value, 
corresponding to the cancellation of Ji, i.e. to the transi- 
tion from AF to FM coupling, between 90 and 95° . Other 
structural parameters can also affect Ji like the angle (3 
between the Cu-0(2) ribbons and the Ge-0(2) segment, 
or the pair distances involved in the exchange path, and 



10 



V. Simonet et al.: Effect of structural distortions on the magnetism of doped spin-Peierls CuGeOs 



more generally, all the parameters that may influence the 
strength of the orbital hybridization. 

5.2 Antiferromagnetism vs distortions 

5.2.1 In-chain magnetic impurities 

Concerning the magnetic impurities on the Cu site, let 
us recall that they drive the system towards an AF LRO 
at a temperature T^r different from that induced by non- 
magnetic Cu-site impurities, slightly lower for Ni and larger 
for Mn. From our estimation of the local a angle (Ta- 
ble III), the magnetic local coupling J|°'^ between Ni^+ 
(spin 1) and its Cu^+ neighbors in Ni-doped CuGeOs is 
expected to be still AF but smaller than the coupling be- 
tween Cu neighbors (Ji = 120 — 1%QK). Both an analysis 
of the thermal dependence of the magnetic susceptibility 
[5] and a resonance model for ESR results ,l3i have in- 
deed evidenced the formation of Ni-Cu pairs with an AF 
coupling J{°^ weaker than Ji (most probable value of J'f'^ 
estimated to be around 75 K [13]). The spin 1 of the Ni im- 
purity and the released spin 1/2 of the Cu therefore form 
an effective spin 1/2 solitonic defect at low temperature. 
In the case of Mn-doping, the a value for the Mn-0(2)- 
Cu path (Table III) suggests that although still AF, 
becomes very weak, which could be reinforced by the frus- 
trating J2 next nearest neighbor interaction. Here again, 
some ESR measurements [8] support this result with the 
observation of an hyperfine spectrum of the ^^Mn (nu- 
clear spin 1=5/2) with fully developed 6 lines for very low 
Mn concentration (0.1 %). An upper value of ~4 mK 
can be estimated from the splitting of the hyperfine lines. 
The detection of quasi free 5/2 spins between Tsp and 
T/v thus demonstrates that the effective coupling between 
Mn^+ and Cu^+ is close to zero. Hence, the onset of the 
AF LRO is mainly due to a collective behaviour of the 
spins 1/2 solitonic defects, decoupled from the free 5/2 
spins of the Mn impurities. It is indeed excluded for the 
direct Mn exchange to be responsible for the transition to 
AF LRO, due to the very large distances between the Mn 
atoms. 

Therefore, the AF LRO transition observed for all the 
in-chain impurities presently studied seems to be mainly 
due to the cooperative behavior of soliton-antisoliton pairs 
(one released per impurity) carrying a spin 1/2, whatever 
the magnetic state of the impurity. Now, another param- 
eter can be considered in order to explain the differences 
between the various T/v(?/) transition lines. The present 
study has revealed differences in the amplitude of the dis- 
tortion around each kind of impurities related to their 
ionic radii. For instance, the first neighbor distances to 
the impurity increases with respect to the Cu environment 
following the sequence Ni/Mg/Zn/Mn (section l4.ip . The 
strongest local distortion, in the case of Mn, clearly affects 
further neighbors (more than two Cu-Cu distances). It is 
interesting to note that the T/v values somehow vary in the 
same way. The distortions induced by the impurities along 
the Cu chains can be very large (Table II) compared to the 
atomic displacements produced by the dimerization itself : 



0.2% relative variation of the first Cu-Cu distance. In con- 
sequence, they will destroy locally the in-chain dimeriza- 
tion and will render the energetic criterion for inter-chain 
phase coherence ineffective. This could in turn affect the 
spatial distribution of the soliton-antisoliton pairs along 
the chain axis, i.e. their overall size or/and the mobility 
of the antisoliton. It would become easier for the antisoli- 
ton to accommodate its distance from the impurity and 
the global size of the soliton-antisoliton pair could be en- 
hanced. As a result, this additional degree of freedom may 
favor the overlap of the AF correlated regions and thus in- 
crease the Tat. This mechanism provides an explanation 
for the reported higher T^r value for Mn with respect to Zn 
due to the larger structural distortion of its environment. 
On the contrary, for Ni, the distortions are the smallest 
implying the lowest T/v, as observed. Furthermore, the 
AF coupling between the spin 1 Ni and the spin 1/2 Cu 
may pin the global spin 1/2 soliton-antisoliton pair closer 
to the impurity. As a result, the overlap between the AF 
correlations is reduced and the T/v is lowered. Lastly, close 
enough to T^r, one can not rule out an indirect polarisa- 
tion of the AF correlated regions by the Mn^+ spin (via 
the inter-chain interactions for instance), which could ac- 
celerate the transition towards the AF LRO. 



5.2.2 Out-of-chain Si impurity 

An unsolved issue remains concerning the strong influence 
of the out-of-chain Si-doping on the magnetic properties of 
CuGeOa. The present study has allowed us to determine 
the large local decrease of the superexchange angle a due 
to Ge substitution by the much smaller Si ion. a was found 
to be smaller than the one calculated in the case of Mn- 
doping for which an almost vanishing was inferred 
(Table III). This strongly suggests that J{°'^^ coupling the 
spins of each Cu pair on both chains facing the Si impurity, 
is indeed ferromagnetic. This is in agreement with theo- 
retical predictions of a weakly ferromagnetic value of J^"*^, 
which was calculated using realistic structural parameters: 
q;=95° and Si-O bond length smaller by 0.13 A than the 
Ge-0 one [TH1[3H]. The authors have invoked additional 
parameters enhancing this trend towards ferromagnetism: 
the further weakening of the hybridization of the 2py 0(2) 
orbitals with Si, with respect to Ge, and the influence of 
frustrating J2 next nearest interaction. This result is also 
comforted by the structure and magnetic properties of an 
isostructural CuSiOa single crystal [23^ 24 . The a angle 
is equal to 94-95° in this compound [23|l40j. i.e. close to 
the one determined locally around Si impurities in doped 
CuGeOa (cf. section IV. B). The analysis, from magne- 
tometry measurements and a neutron powder diffraction 
experiment, of the peculiar AF LRO arising around 8 K 
in CuSiOa, suggests a weak ferromagnetic nearest neigh- 
bor exchange interaction Ji and a stronger antiferromag- 
netic next nearest neighbor interaction J2, similar to the 
CuGeOa one [23,40]. 

Important consequences are expected from a ferromag- 
netic J{°'^ on the nature of the solitonic defects produced 
by the Si impurity and thus on the associated magnetic 
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properties. The dimerization pattern imposes that the Si 
impurity faces one strong bond corresponding to a dimer 
and one weak bond between two dimers. One can rea- 
sonably assume that this local variation of Ji will alter 
only the strong bond, freeing two spins 1/2, and will have 
no effect on the other chain (even locally reinforcing its 
dimerization). In this case, the two released spins 1/2 of 
the dimer, ferromagnetically coupled, may combine into 
an effective free spin 1. 

The Curie constant, extracted from the temperature 
dependence of the susceptibility between Tsp and T^r, in 
a series of CuGci_3;Sia;03 samples [7], was previously an- 
alyzed according to Eq. [3] in a simple model where local 
spins 5=1/2 are associated with the solitonic defects. 



\para (-^) — ^pa 



S{S + 1) 



(3) 



with g = 2.05, the Lande factor, /is the Bohr magneton, 
ks the Boltzman constant, 9 the Curie- Weiss tempera- 
ture, and Kpara = N X with N the number of released 
effective spins S per impurity. Once the contribution of 
spins 1/2 due to intrinsic impurities has been subtracted, 
the number TV of spins 1/2 per Si is found close to 3, a 
puzzling value El- Assuming now that the solitonic defects 
are ferromagnetic pairs of spin 1/2, the susceptibility of 
one pair is given by : 



Xpair- 

(T) 



(.9 Ms)^ 



ksT 3 + eM~Ji/kBT) 



(4) 



At low temperature, the expression of XparaiT) then re- 
duces to Eq. [3] for N such magnetic entities with an effec- 
tive spin 5=1. This analysis yields iV ~ 1.2, which is close 
enough to the expected value of one ferromagnetic pair of 
spins 1/2 per Si impurity in this new description. This can 
be further argued from the analysis of the magnetization 
measured in a 0.3% Si-doped CuGeOa sample (see Fig. 2 
of Ref. 01]). The M{H) curve reaches a first saturation 
plateau Msat = NxgfisS = 0.0077/i_B in the SP phase 
before entering the incommensurate phase. The analysis 
of the whole M{H) curve with Brillouin functions, after 
correction for the intrinsic impurities contribution, yields 
N ~ 1.1 spins 1 or TV - 2.2 spins 1/2. The first result is 
in better agreement with the expected number of effective 
spins 1 and with the previous estimation of N , deduced 
from susceptibility analysis. The revised analysis of the 
magnetic measurements in Si-doped compounds is there- 
fore in favor of the release of one solitonic defect per Si 
with an effective spin 1, on one of both adjacent chains. 

A consequence of the effective spin 1 of the solitonic de- 
fect produced by Si will be to renormalize cx S{S + 1) 
by a factor of 2.67 as compared to the case of spin 1/2 
solitonic defects released by the same concentration of Cu- 
site impurities. However, this comparison only holds if the 
spatial extension of the solitonic defect is otherwise iden- 
tical, thus producing the same overlap between AF cor- 
related regions. The overall size of the soliton-antisoliton 



^ Note that in Ref. [7], the values of Kpara were not cor- 
rected for the 0.11% intrinsic spin 1/2 impurities found in pure 
CuGeOs and expected to be also present in doped samples. 



pair along the c— axis was indeed estimated to be the same 
for Si and Zn from X-ray measurements [42j . This sup- 
ported the origin of the factor ^^3 between the Tjv of Si 
and of non-magnetic Cu-site impurities from the differ- 
ent types of the effective spin associated to the solitonic 
defects, at least for low doping levels (see Fig. [2]). The 
reason why this scaling does not hold at higher x may 
be due to a different mechanism involved in the onset of 
AF LRO. The distances between impurities are becoming 
shorter than the ID magnetic correlations, and moreover, 
a first-order transition is observed in the case of Cu-site 
impurities between a dimerized and a uniform AF LRO, 
whereas the coherence length of the spin-Peierls ordering 
decreases progressively in the Si doping case [i^ . 



5.3 Spin-Peierls transition vs distortions 

Contrary to the sensitivity of the AF LRO to the nature of 
the impurity, discussed in the previous section on the basis 
of the EXAFS results, the decrease of the spin-Peierls tem- 
perature appears to be more universal. On the Cu site, the 
decrease of Tsp has exactly the same concentration depen- 
dence for all studied Cu substituents, whatever the spin 
impurity and the amplitude of the in-chain local distor- 
tions around impurities (see Fig. For Si-doping on the 
Ge site, Tgp decreases ^ 3 times faster with the concen- 
tration, although the size of the solitonic defect induced 
by Si in the chain direction is similar to the one induced 
by Zn impurities [42j . 

Moreover, the average correlation length of the dimer- 
ized regions along c does not seem to be so relevant for 
the appearance of the SP LRO, as shown by the mea- 
surements of fluctuations in pure and doped CuGeOa for 
different impurities above the SP transition [44ll42]. First, 
the regime of ID pretransitional fluctuations along the 
chain direction in pure CuGcOs was shown to start well 
above Tgp, while these fluctuations become 2D and then 
3D much closer to Tsp- Secondly, the average size of the 
dimerized regions along c in the pretransitional regime for 
Si, is comparable to that of Zn despite the factor ^ 3 
in their Tsp- Last, the extension of the ID pretransional 
fluctuations along c is much smaller for Ni than for Zn 
or Mg [4? , whereas their decrease of Tsp has exactly the 
same doping level dependence. 

Therefore, all these experimental results point out that 
the relevant parameter for the decrease of Tsp through 
doping is not the way the dimerization is affected along 
the chain direction, but perpendicular to it through the 
stabilization of the out-of-phase 3D dimerization pattern. 
This is supported by a major difference observed concern- 
ing the spatial extension of the perturbation of the dimer- 
ized regions when doping inside and outside of the Cu 
chains: Only one chain is affected for in-chain doping by 
Zn instead of three chains in the case of out-of-chain Si- 
doping [32] • This has provided Pouget et al. [32] with a 
simple model explaining the factor 3, between Si and Zn 
doping, in the respective decrease of Tgp and in the crit- 
ical concentration above which the SP LRO is destroyed. 
According to the present study, we can generalize this last 
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result concerning Zn to the other Cu substituents, Ni and 
Mg, which produce smaller local distortions. As concerns 
the Mn impurity, which leads to stronger distortions than 
Zn, the fact that Tsp has the same concentration depen- 
dence tends to show that Mn-doping affects only one Cu 
chain. Note that the furthest analyzed atomic distance 
probed by our EXAFS experiment is shorter than that to 
the next Cu on a neighboring chain (4.25 A). In summary, 
the universality of the SP transition seems to be related 
to the transverse weakening of the dimerization which de- 
pends therefore essentially on the position of the impurity, 
outside or inside the chains. 

Lastly, let us come back to the difference in the dis- 
tortions induced by different concentrations of the same 
impurity. The Ni-0(2) distance to the first neighbors in- 
creases more for 6% Ni than for 2% Ni doping with respect 
to the Cu environment. The higher doped compound, also, 
has a more disordered environment. It is interesting to 
note that these concentrations lie before and after the 
end of the Tsp line, and on each side of the transition 
from dimerized to uniform AF LRO. This result under- 
lines some influence of the disorder in the disappearance 
of the spin-Peierls phase. 

6 Conclusion 

In the present study, the EXAFS technique was used as 
a powerful probe of the local atomic arrangement around 
various impurities in the paramagnetic state of doped CuGeO 
This technique has allowed us to quantify the distortions 
in the impurity environment, such as the modification of 
the shape of the O octahedron aroimd the Cu site, and 
to evaluate the local modification of the angles and dis- 
tances relevant for the magnetic properties. An increase 
of the first neighbor distance to the Cu site was observed 
when substituting Cu by Ni, Zn and Mn, in this order. An 
elongation of the distances to the further neighbors was 
evidenced in the case of Mn-doping which thus produces 
the largest alteration among the Cu site impurities. The 
angle a of the Cu-O-Cu superexchangc path leading to the 
AF Ji interaction between Cu^+ first neighbors was found 
to decrease progressively when doping with Ni, then with 
Zn or Mn, and at last with Si. These results have allowed 
establishing the local weakening of the Ji interaction for 
Ni/Cu substitution, its almost cancellation in the case; of 
Mn/Cu substitution and, lastly, its change of sign for Si 
doping on the Ge site, leading to a small ferromagnetic 
coupling between each set of Cu pairs of two adjacent 
chains. 

In a global scheme where each impurity is considered 
to release magnetic solitonic defects, an analysis of the 
structural distortions of the doped CuGeOs in relation 
with their magnetic properties could then be proposed. 
A good description of the magnetic behavior of Si-doped 
compounds is obtained calling upon ferromagnetic pairs 
of spins 1/2, which yields in particular a simple explana- 
tion for the renormalization of Tn for Si with respect to 
Zn-doping. Another interesting consequence of this study 
concerns the Cu-site impurities for which the amplitude 



of the structural distortion, rather than the magnetic mo- 
ment of the impurity, is found to be correlated to the 
variation of Tjy. This last result underlines the importance 
of the structural distortions in the magnetic behavior of 
spin-Peierls doped compounds and opens the way to mod- 
els taking into account realistic structural parameters. 

We would hkc to thank B. Canals, P. Monod, R. Mehn, 
D. Khomskii and N. Laflorencie for fruitful discussions. 
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